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We demonstrate that long-range interaction in a system can lead to a very strong interaction 
between long-wavelength quasiparticles and make them heavily damped. In particular, we discuss 
magnon spectrum using 1/S expansion in 3D Heisenberg ferromagnet (FM) with arbitrary small 
dipolar forces at T <C To- We obtain that a fraction of long- wavelength magnons with energies ek < 
T has anomalously large damping Fk (Fk/ek reaches 0.3 for certain k). Fhis effect is observed both 
in quantum and classical FMs. Remarkably, this result contradicts expectation of the quasiparticle 
fvj ' concept according which a weakly excited state of a many-body system can be represented as a 

collection of weakly interacting elementary excitations. Particular materials are pointed out which 
are suitable for corresponding experiments. 
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I. INTRODUCTION 

The concept of elementary excitations or quasiparticles is one of the most powerful tools in discussion of low-energy 
properties of strongly interacting many-body systems*^ According to this concept a weakly excited state of a system 
can be represented as a collection of propagating weakly interacting quasiparticles which carry quanta of energy e^ 
and momentum k. Elementary excitations, being wave-packets of stationary states, have finite lifetime (or damping 
Tk) which is interpreted as a result of quasiparticles spontaneous decay and interaction between them (at T ^ 0). 
According to the quasiparticle concept supported by quite a general line of argument^— (not considering, however, 
long-range interactions in a system), long- wavelength elementary excitations are well-defined (i.e., £k 3> Tk). It is 
usually the case also for short-wavelength quasiparticles but a limiting number of systems is known in which their 
lifetime is very short or zero. Thus, the spectrum of short-wavelength elementary excitations in liquid He crosses a 
two-particle continuum at threshold momentum k c . As a result spontaneous decay of an elementary excitation into 
two quasiparticles is allowed at k > k c by energy and momentum conservation laws written as^ 

ek = e q + ek-q. (1) 

Decay processes are so strong in liquid 4 He that k = k c is a termination point of the spectrum. A similar behavior 
of short-wavelength quasiparticles (magnons) has been observed recently in a number of magnetic systems: quasi-2D 
spin liquid 5,6 , quasi-lD gapped spin system 6 -^ and in quasi-2D antiferromagnet (AF) with S = 5/2 in strong magnetic 
field 8 . The one-magnon branch disappears completely at k > k c in the quasi- ID material while the ratio r^/ck 
amounts to 0.1 above the threshold in the quasi-2D systems. 

We demonstrate in the present paper that in contrast to the quasiparticle concept expectation even small long-range 
interaction in a system can lead to a very strong interaction between long-wavelength quasiparticles and make them 
badly defined. In particular, we show that a fraction of long-wavelength magnons with ek < T is heavily damped 
(rk/ek reaches 0.3 for certain k) in 3D Heisenberg ferromagnet (FM) on a cubic lattice with arbitrary small dipolar 
forces at small temperature T <c Tq which Hamiltonian has the form 
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where Q^ = {g[J-) 2 (3Ri m Ri m — fippRfnJ/Rfm 1S tne dipolar tensor. Although this model was extensively studied 
before and it describes well a class of magnetic materials (see below) , renormalization of the magnon spectrum by 
thermal and quantum fluctuations has not been analyzed thoroughly yet. It is well known that the spectrum in the 
linear spin-wave approximation (classical spectrum) at L^ 1 C^Cl has the formal 



eok = \J {Dk 2 + g/iH®) {Dk 2 + flf/iJfW + Su sin 2 6» k ), (3) 

uj Q = 4tt(^) 2 , (4) 

where we set the lattice spacing equal to unity, luq is the characteristic dipolar energy, D is the spin-wave stiffness 
which is equal to S J for cubic FM with exchange coupling J between nearest neighbors only, #k is the angle between 
magnetization and k, flw is the intrinsic magnetic field which is zero, e.g., in the multidomain sample 11 and L is the 
characteristic length of a given domain. We assume below that T <C SD ~ Tc (i.e., we do not consider frustration 
which can reduce D considerably and assume that Tc ~ SD ~ S 2 J) and D 3> Sojq as it usually is. 

Long-wavelength magnons are well defined in the model ([2|) with Qfl, — at T <C Tc which damping has the form 
Tk ex T 2 k 4 hik£iH Notice that this previous finding is in full agreement with the quasiparticle concept because 3D 
FM is weakly excited at T <C Tc- We demonstrate below that dipolar forces, despite their smallness and due to their 
long-range nature, give rise to great renormalization of the bare spectrum (J3J. Renormalization of the real part of the 
spectrum has been discussed in our previous paper—. We show there that fluctuations lead to the gap in the spectrum 
which resolve problems of infrared singularities of some observables obtained by other authors (see discussion below) . 
We turn in the present paper to calculation of the magnon damping. 

Dipolar forces lead to three-magnon vertexes giving rise to processes of magnon decay ([T]). However, we obtain 
below that confluence processes reduce the magnon lifetime much stronger. They arise only at T ^ and have the 
following conservation law (cf. Eq. (pj): 

£k = e q - £k- q - (5) 

The important role of these processes for magnon relaxation was recognized long time ago in Refs. 14 . These works were 
motivated by non-linear ferromagnetic resonance experiments. Therefore large external magnetic field was taken into 



account there leading to g/iH^ > Sloq- It is seen from expressions obtained in RefsJ^ that £k 3> Tk at g/iH^ ~ Sluq 
but Tk — > cm as k, H^ — > and results should be reconsidered if H^' — 0. We perform this reconsideration below. 

We derive magnon damping using 1/5 expansion. It is shown that in accordance with previous results^ T^ 
diverges as k — » in the first order in 1/5 that is a consequence of the Goldstone character of the bare spectrum 
([3]) at H^> = 0. This divergence signifies that the spectrum cannot be found using the conventional 1/5 expansion. 
One notice, however, that this singularity should be screened by the gap in the spectrum. That is why the easiest 
way to find the spectrum in this case is to perform the self-consistent calculations. As a result we show that the 
main corrections to the spectrum comes from diagrams of the first order in 1/5 whereas those from higher order 
ones are small by the parameter Sujq/D <C 1. We demonstrate below that magnons are well defined at T = 0, while 
thermal fluctuations lead to a great enhancement of the damping: a peak arises in the ratio I\/ek at small k which 
height reaches 0.3 for momenta directed nearly along magnetization. The fraction of overdamped magnons is small 
and one could expect small influence from them on the system properties (magnetization, specific heat, etc.). By 
reducing dipolar forces radius of action we demonstrate that it is their long-range nature that is responsible for such 
a remarkable suppression of the long-wavelength quasiparticles. We show that this suppression can be seen both in 
quantum and classical FMs because thermal fluctuations are responsible for it. 

Interestingly, quantum and thermal fluctuations lead to smaller I\/ek in lower dimension 2D FM on square lattice 
with dipolar forces that is discussed in our previous paper Ref^.. We find there that thermal fluctuations also lead 
to a peak in I\/ek at small k in 2D FM which height, however, is of the order of T/Tc <§C 1 for 5 ~ 1 and reaches the 
value of 0.16 for S> 1. We show below that the origin of the greater role which play fluctuations in higher dimension 
is that dipolar forces give rise to larger three-particle vertexes in 3D FM. 

The rest of the present paper is organized as follows. The Hamiltonian transformation and the technique are 
discussed in Sec. HU Renormalization of the energy and the real part of the spectrum are discussed briefly in Sec. IIIII 
Magnon damping is considered in Sec. IIV1 In Sec. [V] we i) show that the observed anomaly in the damping can be 
seen both in quantum and classical FMs, ii) compare the damping in 2D and 3D FMs and make a counter intuitive 
conclusion that it is smaller (compared to the real part of the spectrum) in lower dimensional FM, and iii) discuss 
particular compounds suitable for corresponding experimental verification of our results in 3D FM. Sec. IVII contains 
our conclusion. An appendix is included with some details of calculations. 



II. HAMILTONIAN TRANSFORMATION 

After the Fourier transformation Hamiltonian ([2|) is written as 

n = -\ J2 ( J k<W + Qf) szst k - g^Hms*, (6) 

k 

where Jk = X)z Jim exp(zkR; m ) and Q^ = ^2 t Qfm exp(zkR; m ). Notice that we take into account the Zeeman term 
—g\xH 'Y^i Sf in Eq. ([5]) that is necessary to do in order to describe the real finite size samples: H is the field in a given 
domain produced by all other domains in the multidomain sample in zero external magnetic field. We will assume 
that if = in the unidomain sample that is magnetized in the direction in which it can be considered to be infinite. 
The dipolar tensor Q^r possesses the well-known properties ^ 10 ' 16 

Q k -u Q \— — l,if-<fc<l, (7) 

Qf=u ( l-Af a )s a p, (8) 
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where A/" Q are demagnetizing factors. After Dyson-Maleev transformation 
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Hamiltonian © has the form V. = E + J2i=i ^»! where E is the classical ground state energy and Hi denote terms 
containing products of i operators a and a' . In particular, Hi = because it contains only Qo with a ^ j3 and 
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where 9t is the number of spins in the lattice, we drop index k in Eqs. (jlll) and (|12[) and 



£ k = 5(J -Jk)-|(^ k ' 
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where expressions after w are approximate values of corresponding quantities at L^ 1 <C k <C 1 and </>k is the 
azimuthal angle of k. The expression in the brackets of the last term in Eq. (fT3|) is the intrinsic magnetic field H^K 
In the multidomain sample the term 4ngfj,SAf z is the demagnetizing field of the considered domain that is equal to 
H in the domain volume so that H"> = OJ^- The intrinsic field is zero also in a unidomain sample that is infinite in 
the direction of magnetization if the external field is zero because Af z = in this case. One leads to Eq. ([3]) from 
Eqs. (|10[) . (|13|) and (|14p for the spectrum given in the linear spin- wave approximation by 
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To perform the calculations it is convenient to introduce the following retarded Green's functions: G(w,k) = 

(ak,a k )w, F(oj,k) = (a k ,a_ k ) w , G(w,k) = (al k ,a_k) w = G*(-w,-k) and Ft(w,k) = (ol k , 0^)01 = F*(-w,-k). 
We have two sets of Dyson equations for them. One of these sets has the form: 

G(u, k) = G<°> («, k) + G<® («, k)E(«, k)G(w, k) + G<°) (u>, k) [B* + n(w, k)]Ft («, k), 
Ft( w ,k)=G (0) (w,k)E(w,k) J F 1 t( Wjk ) + G (0) ( w ,k)[B k + nt( w ,k)]G(w,k), 

where G^ (w, k) = (u> — E^ + i5)~ 1 is the bare Green's function and E, E, II and II ^ are the self-energy parts. Solving 
Eqs. (fi"6]) one obtains: 

+ £ k + E(w,k) 
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G(w,k) = - 
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where 



P( W ,k) = ( W + *5) 2 _-e 2 k -^( W ,k), _ _ (18) 

fi(w, k) = £ k (E + E) - B k n - B k n f - (w + i<J)(E - E) - IilT + EE (19) 

and eok is given by Eq. (I15|) . The quantity f2(w,k) describing the spin- wave spectrum renormalization is considered 
below. The last two terms in Eq. (|19l) give corrections of at least second order in 1/5. 




FIG. 1. Diagrams for self-energy parts of the first order in 1/S. 

III. RENORMALIZATION OF THE ENERGY AND THE REAL PART OF THE SPECTRUM 

The classical ground state of the model ([2]) is continuously degenerate: magnetization can be oriented in any 
direction. It is well known, however, that quantum fluctuations lead to anisotropic corrections to the energy selecting 
a limiting number of states . ' 18 Intrinsic anisotropy of the dipolar interaction is the origin of this effect. Anisotropic 
part of the first l/S'-correction to the energy E$ has the form 



^C^^hl+lht+llll), (20) 

^V ^ ' (21) 

where 7^ are direction cosines of the magnetization and components of the dipolar tensor in Eq. (|21[) are taken relative 
to cubic axes. The constant C can be calculated numerically and one obtains in accordance with Refsi 17 ' 18 that it is 
positive for simple (C sw 0.012) and negative for face-centered (C « —0.005) and body-centered (C ps —0.04) cubic 
lattices. Then, an edge of the cube is easy direction in the simple cubic lattice whereas a body diagonal of the cube 
is easy direction in face- and body-centered cubic lattices. 

It is well known that fluctuations leading to anisotropic corrections to the energy which low the energy symmetry 
to a discrete one naturally lead also to a gap in the spectrum. To mention only a few, examples are antiferromagnet 
containing two coupled antiferromagnetic sublattices^ and square planar rotator model with dipolar interaction^ 
As we show in Refs i 13 i 15 , it is also the case in the considered 3D FM and in 2D FM with dipolar forces. 

Renormalization of eok in the first order in 1/S comes from two diagrams shown in Fig.[TJ Both diagrams contribute 
to renormalization of the real part of the spectrum and only the loop diagram (b) leads to magnon damping. Renor- 
malization of the real part of the spectrum was discussed in detail in our previous paper Ref J^. It is demonstrated 
there that renormalized spectrum has the form 
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^=Je 2 u + 2—k 2 + A 2 sm 2 k , (22) 



Y ^-> T < y/SDuo, 



A ~< 1 S^IT ^ 



2 8 DyfSD^'' 



T » y/SDu , 



where Y = C/2 f=a 0.006 for simple cubic lattice and Y = — 3C/4 for face- and body-centered cubic lattices. Notice 
that we present in Eq. (|22|) only corrections which change drastically the bare spectrum eok at T <C Tc- The rest 
corrections observed in Ref>^ lead to only small renormalization of constants D and ujq in Eq. (|3|). In particular, the 
term in Eq. (l22l proportional to k 2 makes the spectrum linear at sufficiently small k and sin6>k = whereas eok oc k 2 . 
The last term under the square root in Eq. (|22|) is the square of the gap value. Observation of the gap induced 
by dipolar forces resolves problems of infrared singularities of longitudinal spin susceptibility and corrections to the 
spin- wave stiffness in 3D FM 13 . Eq. (|2"2"|) can be simplified as follows in three limiting cases which are considered 
below: 
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where we assume in the second and the third lines that sin 0k is not too small: sin 0k 3> k^J D / Sojq <C 1. 



IV. MAGNON DAMPING 



Imaginary part of the loop diagram shown in Fig. HJb) can be obtained quite straightforwardly as it was done, for 
instance, in our previous paper— devoted to 2D FM with dipolar interaction. Some details of these simple but tedious 
calculations together with general expression for Im£7(w, k) in the first order in 1/5 can be found in Appendix. The 
value ImJ7(w, k) is an odd function of u> and we calculate it for u> — ek only. According to Eq. (| 18 j) the spin- wave 
damping Tk at momentum k is given in the first order in 1/5 by the relation 



Inrfi(w = £k, k) 
2^ ' 



(25) 



Strictly speaking, one has to put the bare spectrum eok in Eq. (1251) instead of ek in the first order in 1/5. It is shown 
in Appendix that ImJ7(w = eok, k) = at sin 9k = 0. But, as we show below, the damping diverges as k — ¥ at T ^ 
and sin #k 7^ 0. This divergence signifies that the spectrum cannot be found in conventional way as a series of integer 
powers of 1/5. However one may try to find a solution of the equation 2?(u;,k) = self-consistently using general 
expressions for self-energy parts obtained within the first few orders in 1/5. It was the way in which the real part of 
the spectrum ([2"2"|) was obtained in RefJ^. The only assumption we made there is that Tk ^C £k at fc ^ A/y/SDuJo 
(that is really the case as results show presented below). Notice, in particular, that the spectrum (|2"2")l is of the order 
of VS at fc = 0, sin 0k 7^ and T <C y/SDujQ, while the bare spectrum is of the order of 5 1 and corrections to it of 
the first order in 1/5 are of the order of 5°. Such dependence of the spectrum on fractional powers of 5 is the result 
of the self-consistent procedure being used. 

We obtain below Tk self-consistently using the expression for self-energy parts in the first order in 1/5 and assuming 
that £k is the renormalized spectrum given by Eq. (|22[) . In particular, the gap in the spectrum screens the divergence 
of Tk obtained in the first order of the conventional 1/5 expansion. The results presented below for k 3> \J Sujq/ ' D 
are valid for all #k whereas the damping was found at k <C \J ' Sujq/ ' D under assumption that sin #k 3> k \J D j Sujq <C 1 . 
The analysis becomes more complicated at k <C ^Sluq/ D in the narrow interval of angles given by inequality sin #k ^ 
k^D/Suo <C 1. The damping is expected to be small at such 6*k because Tk = at sin 6k = in the first order in 
1/5, as it is noted above. It is shown below that expressions for self-energy parts of higher orders in 1/5 involved 
in the self-consistent calculations give negligibly small contribution to the result by the parameter Sujq/D <1. It is 
convenient to consider two regimes, T — and T 3> 5wo, in which only quantum and thermal fluctuations determine 
the damping, respectively. 



Decay processes ([T]) lead to magnon damping at T — which has the form 
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where £k is given by Eq. (|2"2"j) and it is taken into account that e k ~ Dk 2 at k >• ^Sljq/D. It is seen from Eqs. (|24|) 
and (|26p that Tk is much smaller than ek by the parameter Sujq/ D <C 1 and magnons are well-defined quasiparticles 
at T — 0. It is seen from Eqs. (|24|) and (|26|) that for a given 9u the damping decreases monotonically as k decreases 



in the interval A/^/SDujq <C k -C 1 and it is flat at k <€. A/^/SDluq. 



B. T > Sujq 



As we note above, the damping diverges at T ^ as Tk-j-o ^ Tuj q /eok but one obtains a finite Tk as a result 
of self-consistent calculation. We obtain assuming that T ^> ek after tedious calculations some detail of which are 
presented in Appendix 




FIG. 2. Function /(#k) is shown that is given by Eq. ((28 
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and H(x) is the Heaviside step function. As is seen from its graphic shown in Fig. [2] values of the function /(0k) 
lie between 1.2 and 3.1. Both decay (HJ and confluence ([5]) processes contribute to the damping at k ^> ^/SlJqJD 
whereas confluence and decay processes dominate at JC <C fc <C ^Sluq/D and fc <C /C, respectively. 

It is seen from Eqs. (|23[) . (f24| and (f27| that the damping rises as Tu:f i /(D 2 k) upon the momentum decreasing up 
to k ~ A/\/SDujo, it decreases linearly as \/Dujok at /C <C fc <C A/^/SDuq and it is flat at fc -C /C. Thus there is a 
peak in the damping at fc ~ A/^/SDwq. We draw Tk schematically in Fig. [3] using Eq. (|27|) in the particular case of 
T 3> \/SDojo when A is given by the second line in Eq. (|23|). 



The ratio Tk/fik follows qualitatively the behavior of Tk and has a maximum at k ~ A/y/SDuo too. The peak 
height can be estimated from the third line in Eq. (|27|) at fe ~ A/^/SDuiq that gives Tk/ek ~ T/^/SDuo < 1 at 
T <C \/SDujq and Tk/ek ~ 1 at T > ^/SDujo. The peak height cannot be calculated analytically. We perform 
numerical integration to find Tk/^k at k ~ A/t/SDcJq in the particularly interesting case of T ^> s/SDujq which 
results are shown in Fig. [4] (see Appendix for detail). Asymptotics obtained from Eqs. (J22I) and (|27|) are also shown 
in Fig. |U by dashed lines. As the asymptotic for A/^/SJJlJq « K< -JSwq/D was found with the "logarithmic" 
precision, we have introduced a factor of the order of unity under the logarithm in order to improve the quantitative 
agreement in the near vicinity of the peak between results of the analytical consideration and numerical integration. 
In particular, the factor introduced was equal to 1.03, 1.45 and 2.25 for 0k — tt/2, tt/A and tt/8, respectively. It is 
seen from Fig.|3]that the peak becomes sharper and higher upon decreasing of sin 0k and its position moves to smaller 
momenta. Characteristics of the peak, its height and position, as functions of 0k are presented in Fig. [5] It is seen 
that the value of the peak height increases from 0.185 at 9k — 7t/2 to 0.296 at sin 8k <C 1. The rapid decrease to zero 
at sin 9k < k^D/SujQ <C 1 discussed above is not shown in Fig. [5j 




FIG. 3. Sketch of the damping F^ at T 3> \/jDujq given by Eq. (|27[) . Thermal fluctuations are responsible for the peak at 
k ~ A/\/SDujq: damping at T = given by Eq. (|26p rises monotonically as fc increases. 
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FIG. 4. (Color online). The ratio is shown of the magnon damping Tk and the real part of the spectrum ek at T > \/SDujq 
as a function of reduced momentum k = ky/SDujoj A in the vicinity of the peak for three particular values of #k- Asymptotics 
at k <C A/^/SDujo and k 3> A/\/SDojo are shown by dashed lines which are obtained from Eqs. (|24p and (|27|l . 



DISCUSSION 



Classical-spin limit. — It should be noted that anomalously large damping obtained above can be seen both in 
quantum and classical FMs because thermal fluctuations are responsible for it. Really, we conclude from Eqs. (|22|) 
and (|2"7|l taking into account that Tc oc S 2 J at large S that at fixed ratio T/Tc one has £k, Tk ~ S. Thus both Tk/fik 
and the peak position are S'-independent. As it is explained in Refs . 15 ' 21 ' 22 , the spectrum and corrections to it in the 
classical limit can be obtained from expressions found above by multiplying them by S and taking the limit 
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FIG. 5. (Color online). The peak height (Tk/eijmax and position fc max = K max A / V SDljq of the ratio r^/ek at T> y/SDuio 
as functions of 6k- 



FIG. 6. Some diagrams of higher order in 1/5. 



Quantum corrections (i.e., T-independent ones) die out as a result of this procedure. 

Role of the long-range character of the dipolar forces. — It should be stressed that the long-range nature of dipolar 
forces is the origin of such a large damping obtained. It can easily be shown by restricting dipolar forces radius of 
action to a few lattice spacings. Such a restriction leads to small damping because the value of three-magnon vertex 
dill) becomes of the order of wo/c 2 if fci,2,3 ~^<1 rather than ujq. It is well known that small short-range interactions 
of other types lead to small damping as well. 23 

Comparison with 2D FM with dipolar forces. — Interestingly, quantum and thermal fluctuations lead to smaller Tk 
(compared to ek) in lower dimension 2D FM on square lattice with dipolar forces. It is demonstrated in Ref J£ that 
thermal fluctuations lead also to a peak in Tk/ek at small k in 2D FM which height, however, is of the order of 
T/Tq <C 1 for S ~ 1 and reaches the value of 0.16 for S>1. The origin of the greater role which play fluctuations in 
higher dimension is that Q? x u ok in 2D FM whereas Q? K w o in 3D FM that leads to larger three-particle vertex 
in 3D FM. 

Higher- order diagrams. — The interesting property of the results discussed here is that higher-order diagrams some 
of which are presented in Fig. |6j being included in the self-consistent calculations, give a negligibly small contribution 
by the small parameter Sujq/D. The origin of this fact is that singularities of higher-order diagrams are screened by 
the gap and all vertexes are small being by the order of magnitude not larger than wo when corresponding momenta 
are much smaller than ^JSujq/D. 

Experimental verification. — The main obstacle for experimental observation of anomalous damping in FMs discussed 
above is magnetocrystalline anisotropy which contribution to the gap in the spectrum A a is much larger than A at 
T <C Tq in the majority of materials. The isotropic gap A , in contrast to anisotropic one of the dipolar origin 
A sin #k i leads to exponential decay of the damping 



Tk oc exp 



ATDk 2 



at 



k< 



A„ 



(31) 



screening the effect we discuss. However, there is a class of insulating FM materials which magnetism is due to 
magnetic ions in S-state and which, as a consequence, have very small magnetocrystalline anisotropy. The most 
suitable compounds for corresponding experiments seem EuS (Tc = 16.6 K, S = 7/2, D sa 1.68 K, u; ~ 0.58 K, 
A a k 0.005 Kp. and CdCr 2 Se 4 (T c - 130 K, S = 3/2, D « 3.9 K, u « 0.3 K, A Q w 0.002 K) 25 . It should be 
noted that because Suo > D, asymptotics (|26p and ([2T]) for Tk are not valid for EuS. However, calculations of Tk at 
k ~ A/\JSDujo and sin 0k ~ 1 showing the peak remain valid for EuS. Particular estinrations demonstrate that the 



maximum of the damping in EuS and CdCr 2 Se4 at #k 
respectively. 



n/2 and T = 0.1 -f- 0.2T C is at k ss O.OlA and 0.002A~ 
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VI. CONCLUSION 



To conclude, we demonstrate by the example of 3D Heisenberg FM with arbitrary small dipolar forces at T <C Tq 
that long-range interaction in a system can lead to a very strong interaction between long-wavelength quasiparticles 
and make them heavily damped. We show that magnons are well-defined quasiparticles in 3D FM at T = which 
damping is given by Eq. (|26|) . Thermal fluctuations give the main contribution to the damping at T> Sloq which is 
given by Eq. (|27]l and which is also sketched in Fig. [3] There is a peak in the damping at such T at k max ~ A/^/SDujq. 
There is also a peak at k ~ k max in the ratio T^/e^ which height becomes of the order of unity when T 3> \/SDuq. 
The ratio Tk/ek is shown in Fig. 3] in the vicinity of the peak. The peak height and position are plotted in Fig. [5) The 
maximum value of Tk/ek approaches 0.3 at small sin#k- It should be noted that only small fraction of magnons with 
£k < T appears to be overdamped. That is why one could expect small influence from them on the system properties. 
However, the observed suppression of long-wavelength magnons is remarkable because it contradicts expectation of 
the quasiparticle concept. 



Because thermal fluctuations are responsible for the peaks in Tk and Tk/ek) they can be seen both in quantum and 
classical FMs. We demonstrate that it is the long-range nature of dipolar forces that is responsible for the anomalously 
large damping observed in 3D FM: magnons are well-defined quasiparticles if one restricts dipolar forces radius of 
action to a few lattice spacings. We make a counter intuitive conclusion that dipolar forces lead to smaller magnon 
interaction in lower dimension 2D FM discussed in our previous paper—, where the peak in the ratio Fk/ek is also 
observed which height, however, is of the order of T/Tc <C 1 for S ~ 1 and it reaches 0.16 for S 3> 1. We show that 
magnetocrystalline anisotropy leading to isotopic contribution to the gap A a screens the enhancement of the damping 
we discuss because, in contrast to anisotropic gap A sin #k of the dipolar origin, it leads to exponential decay of the 
damping at small k (see Eq. (|31[0 . Seemingly surprising fact that such a remarkable anomaly of the damping was not 
obtained before can be explained by small value of k max and by large value of A a as compared to A in the majority 
of materials. We show that this effect can be observed experimentally in materials with small magnetocrystalline 
anisotropy the most suitable of which seem EuS and CdCr2Se4. 
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Appendix A: Details of the damping calculation 

Taking into account that summation over small momenta is essential in calculation of the damping at k <C 1 and 
T <C Tq one has for the contribution to Q(iu>, k) from the loop diagram shown in Fig. [TJb) 

x {sin 2 29 ii (B 1 B; +E X E 2 + Wicj 2 ) (Ala) 

— 2 sin26>k(Sw -Dq 2 sin2#i sin 2 # 2 sin(0i — 2 ) sin(0 2 — 0k) 

- (Dq 2 Dq 2 + u>iw 2 ) sin26q cos(0i - k )) (Alb) 
+ (2Dq 2 Dq 2 sin 2 20 x + Suj Dql sin 2 2 sin 2 26>i 

+ {Dq\Dql - wiu; 2 ) sin2#i sin26> 2 cos(0i - 2 ))} (Ale) 



■sm 2 k T JJ £, [(iu x ) 2 - e 2 ][{iuj 2 ) 2 ~ e 2 ] 

x {Dq\Dq\ sin 2 20i + Su Dq\ sin 2 2 sin 2 20i cos 2 (0 2 - k ) 

+ Dq\Dq 2 sin20i sin20 2 sin(0i — 0k) sin(0 2 — k ) 

— ujiuj 2 sin20isin20 2 cos(0i — 0k)cos(0 2 — 0k)} (Aid) 



zw- 






29T * f- f- [(* Wl ) 2 - e 2 ] [(i W2 ) 2 - e 2 ] 



x {£>q 2 sin20 2 (sin20 2 + sin20i cos(0i - 2 )) (Ale) 

+ sin2#k(-D<Z 2 sin20 2 cos(0 2 — 0k) — Dq 2 sin20i cos(0i — 0k) 

— Sujo sin 2^i sin 2 2 cos(0 2 — 0k) cos(0 2 — 0i))}- (Alf) 

We obtain after summation over imaginary frequencies and analytical continuation on to from an imaginary axis to 
the real one 



\ "1 



1 



[(iuix) 2 - e 2 ][(iu 1 - iuj) 2 - e 2 } 



-?— (1 + Ai + JV 2 )sgn(u;)<$(|u;| - ei - e 2 ) - — ^— (Ai - N 2 )(S(uj - e x + e 2 ) - 6{lu + ei - e 2 )), (A2a) 

4eie 2 4eie 2 

(iwi)(iw — iwi) \ 



Im r£ 



[(zwi) 2 — e 2 ][(iwi — iuj) 



^(1 + Ai + A 2 )sgn(w)5(M - e x - e 2 ) + ~(Ai - A 2 )(<5(w - e x + e 2 ) - S{uj + ei - e 2 )), (A2b) 



Im T E 



ZWi 



[(iwi) 2 — e 2 ][(iwi — iw) 2 — 



-?-(l + Ni + N 2 )6(\oj\ - ei - e 2 ) - -^-(Ai - A 2 )(5(w - e x + e 2 ) + 5{u + e x - e 2 )). (A2c) 

4e 2 4e 2 



The first and the second terms in Eqs. (IA2[) describe magnon decay and confluence processes, respectively. It can be 
straightforwardly shown using Eqs. (|A1|) and (|A2|) that Im£7(u;,k) = at sin 0k = if e lt2 and ek are bare spectra 
given by Eq. ([3]). 

1. T = 

Only decay processes contribute to the damping at T = because Ai = A 2 = in Eqs. (JA2I) . The main contribution 
at 1 > fc > ^/Sluq/D comes from terms (c) and (e) in Eq. (|A1[) . Summation over momentum q\ ~ fc is essential in 
this case that leads to the first line in Eq. (f26| . 



The main contribution at A; ^; y/Suo/D comes only from term (d) in Eq. (|A1|) . Summation over q\ ^> k and 
sin 6>i <C qiy/D/(Suj ) is essential that leads to the second line in Eq. (|2^|) . where £k should be replaced by eok in the 
first order in 1/S. 
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Although there is no problem with infrared singularities at T = 0, one has to go beyond the first order in 1/S to 
obtain correct expressions for the damping. It is the consequence of the fact that the bare spectrum is renormalized 
greatly by the first 1/S corrections at small momenta. The easiest way to do it in the present case is to find the 
damping self-consistently by using "dressed" Green's functions in the loop diagram and showing that higher order 
diagrams are negligible some of which are shown in Fig. [6j Self-energy parts in numerators of Green's function (|17p 
are negligible at k -C 1 and T <C Tc and one can use Eq. (|A1|) for the self-consistent calculation assuming that ei,2 
and ek are renormalized spectra given by Eq. (|22|) . As a result one leads to Eq. ([26]) for the damping. 



2. T > Suo 

To simplify the consideration we assume that T 3> e^ that allows to replace N\^ in Eqs. (|A2[) by T/t\fl to obtain 
the main contribution to the damping. Terms (c) and (e) in Eq. (|A1[) play the main role at k 3> ^/Suig/D and one 
leads to the first line in Eq. (|2"T|) after self-consistent calculations. Notice that Im£l(u;, k) has a logarithmic singularity 
in the first order in 1/S (the factor \n(SuJo/A) diverges in Eq. (|2"T|) as A — >• 0) that is screened by the gap in the 
self-consistent procedure. 

The main contribution at k <C y^Stoo/D comes from term (d) in Eq. (|Alj) . It can be readily shown that the energy 
conservation law of the confluence process ([5]) cannot be fulfilled in the limit of small k due to the gap in the spectrum. 
It is easy to realize that the corresponding contribution becomes exponentially small at k < A/\/TD. In contrast, 
equality describing magnon decay ([1} is fulfilled by momenta which are nearly parallel to magnetization: ij> k and 
sin0 q -c qy / D/(Sujo). Confluence processes make the main contribution to imaginary part of term (d) in Eq. (|Alj) 
at A/^/SDujq <fcc y/ScJo/D. But their contribution decreases at k < A/y/SIkJo as k decreases and it becomes of 
the same order as that from the decay processes at k ~ /C = max{A 5 / 2 /(S 2 ujQy/~D), A/ y/TD}. Term (d) in Eq. (|Al|) 
can be represented in the form at K, <C k <C tJSujq/D 

4A^ fl f M f°° A sin 2 2fl(2q 2 + ,4 2 (sin 2 fl-2 g 2 )) 

7 {2q 2 + l)(q 2 + sin 2 6») 3 / 2 Vl - A 2 
\/l + k 2 \Jq 2 + sin 2 9 2q 2 cos9 

(2q 2 + l)ain6 + {2q 2 + 1) sin 



r k = 4_ s in0 k / d9 / dg — ^* ^ , j ; H " H(l-A 2 ), (A3) 

' - ' * '■" 2 + l)(g 2 +sin 2 r° ' ■' " 

V i •- n y '/ — : ' > "> " 2o cost) . 

A = -. ^-TT^^ + CT TT^T^ cot^ (A4) 



where k — ky/SDuo/A and H(x) is the Heaviside step function. Integrals in Eq. (|A3|) can be taken if k ^> A/y/SDuo 
and fc <C A/y/SDuiQ that leads to the second and the third lines in Eq. (|27l) . respectively. The peak height which 
is located at k ~ A/y/SDuo cannot be calculated using Eq. (|A3[) analytically. Results of numerical integration at 
T > y/SDu! using Eq. (|A"3"|) are shown in Fig. g) 

Decay processes come into play at k <C K in term (d) of Eq. (|A1|) , where summation over gi 3> fc and sin B\ <C 
qi \J D / {Sloq) is essential. One leads to the last line in Eq. (f27|) that is simply the last line in Eq. ([26]) multiplied by 
4T/e k . 
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